Mode-locking transitions in nano-structured weakly disordered lasers 
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We report on a statistical approach to mode-locking transitions of nano-structured laser cavities 
characterized by an enhanced density of states. We show that the equations for the interacting 
modes can be mapped onto a statistical model exhibiting a first order thermodynamic transition, 
with the average mode-energy playing the role of inverse temperature. The transition corresponds 
to a phase-locking of modes. Extended modes lead to a mean-field like model, while in presence 
of localized modes, as due to a small disorder, the model has short range interactions. We show 
that simple scaling arguments lead to observable differences between transitions involving extended 
modes and those involving localized modes. We also show that the dynamics of the light modes 
can be exactly solved, predicting a jump in the relaxation time of the coherence functions at the 
transition. Finally, we link the thermodynamic transition to a topological singularity of the phase 
space, as previously reported for similar models. 
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I. INTRODUCTION 

Laser mode-locking (ML) is well known in standard op- 
tical resonators, which are characterized by equi-spaced 
resonances p], Q ■ ML in such a kind of systems is a valu- 
able route for the generation of ultra-short pulses, in par- 
ticular when it is "self-starting" , as due to the nonlinear 
interaction between laser longitudinal modes [3|- Given 
the growing interest in high-Q microresonators and pho- 
tonic crystals 0, @, it is interesting to consider ML in 
integrated devices, which could trigger a new generation 
of highly-miniaturized lasers emitting ultra-short pulses 
(see e.g. Q). 

In this respect, there is a remarkable difference be- 
tween standard resonators and nano-structured cavities; 
indeed the latter are characterized by a non-uniform dis- 
tribution of resonances, given by a strongly modulated 
density of states (DOS) [5J. This situation favors a new 
formulation of the analysis of the self-mode locking tran- 
sition, based on a mean-field thermodynamic approach: 
this is the topic of the present paper, also including the ef- 
fect of some disorder in the system. The thermodynamic 
approach to multi-mode interactions in various physical 
frameworks is well established Q- For example, it was 
recently applied to transverse-mode interaction in res- 
onators [1,0, as well as to standard-laser mode- locking 
transition [nj EH E3] • This transition can be described 
in terms of an effective temperature T, which encom- 
passes the level of noise due to spontaneous emission and 
the amount of energy stored into each mode. At high 
T the mode-phases are independent and rapidly varying 
("free-run" or "paramagnetic phase"); conversely, either 
reducing the spontaneous emission noise or increasing 
the pumping rate, a low-temperature ("ferromagnetic") 
phase can be reached, corresponding to the mode-phases 
locked at the same value. 

A paradigm that has been recently shown to be very 
effective for describing the nonlinear interaction of many 



"modes", and the resulting phase transitions and/or ki- 
netic arrest is the potential energy landscape (PEL) ap- 
proach (see e.g. [13, [I3])- The PEL, as a manifold in the 
configurational phase space, has many stationary points 
(typically minima and saddles) [To] ], whose distribution 
strongly affects the thermodynamics (and the dynam- 
ics) of the system. Recently this paradigm, developed 
to investigate the glass transition phenomena, has been 
applied to the field of photonics, including optical soli- 
tons [H, [TtJ and random-lasers [H, Ojj]. In this respect, 
it is worth to note that the geometrical interpretation 
of the laser threshold was recognized since the begin- 
ning of laser theory, and is considered as one of the suc- 
cessful applications of catastrophe-theory, which classifies 
the singularities of multi-dimensional manifolds @, [20| . 
It is not surprising, therefore, that the mode locking 
transition can be interpreted according to the thermo- 
dynamic/topological transition point of view. Extending 
the topological approach to the nano-lascr is interesting 
for different reasons: on one hand this provides an ele- 
gant and comprehensive theoretical framework to laser 
theory, on the other hand it can be relevant from a fun- 
damental physical perspective. Indeed, in recent litera- 
ture the link between geometr y an d thermodynamics has 
been strongly debated [H [H, 1H, HI , and it is still to be 
established if this theoretical circumstance has physical 
consequences. Thus it is important to identify physical 
systems which can be treated by analytical solvable mod- 
els (for what concerns both the thermodynamical and the 
topological properties) . 

Here we show that mode-locked laser nano-cavities fall 
within this category. In addition we also report on ex- 
plicit expressions of the first-order coherence function of 
the laser emission. This analysis predicts a jump in the 
relaxation time (and correspondingly in the laser line- 
width) at the mode-locking transition. The scaling prop- 
erties of the threshold average mode-energy at the tran- 
sition are found to be strongly sensible to the degree of 
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localization of the involved modes. 

The outline of this manuscript is as follows: in section 
Hi! we will recall the mode-coupling approach to multi- 
mode lasers; in section IIIII we will discuss the physical 
signatures of transitions involving either localized or de- 
localized modes; in section llVl we will report on the ther- 
modynamic approach; the analysis of the topological ori- 
gin of the laser transition is given in [Vj in section I VII we 
briefly discuss the solution of the dynamics of the model; 
conclusions are drawn in section fVIII 



II. MULTI-MODE LASER EQUATIONS 

The coupled mode theory equations in a nonlinear cav- 
ity can be written in the form [ll], [2^, [2(| 

^IT = ~~T~T _ a s a s (t) + (j s - g s \a s {t)\ 2 )a s {t) + r) s (t) 
at oa* 



dU 



with 



and 



H = H + H! 



(1) 



(2) 



S 



g s \a s 



(3) 



In s = 1, 2, .., N with N the total number of modes, 
while a s is the complex amplitude of the mode at oj s , such 
that £ s = u> s \a s \ 2 is the energy stored in the mode. Radi- 
ation losses and material absorption, are represented by 
the coefficient a s , while j s — g s |a s | 2 represents the sat- 
urable gain term and, as usual, the quantum noise term 
due to spontaneous emission is given by a random term 
rj s such that (rj s (t)ri p (t')) = 2k B T bath 8 sp 8(t - t'), where 
ks is the Boltzmann constant and Tb a th is an effective 
temperature (l9l ] . The nonlinear interaction term is 



H, = Rel 



■1 r 

. / j gspqrO, s O,pa q a r \, 
{spqr} 



(4) 



where the sum is extended over all mode resonances such 
that uj s + Lo p = uj q + uj r . The term s = p = q = r is not 
included as it is already described by g s ; the Hamiltonian 
Hj describes mode interaction due to the nonlinearity of 
the gain medium. The field overlap is given by 



^/u s ujpUj q uj r r 

9 spqr = ^ J Xa/3~tS(U s ; OJq, U r , — LO p , T)X 

x EfE^E^E s r dV . 



(5) 



where V is the cavity volume, x 1S t ne third order- 
susceptibility tensor due to the resonant medium and E p 
are the components (a = 1,2,3) of the vectorial mode 
of the cavity at the resonance u) p . \ is give n, in the 



simplest formulation, by the Lamb theory [25j, |26( and, 
neglecting mechanisms like self- and cross-phase modu- 
lation (which give phase-independent contribution to the 
relevant Hamiltonian, see below) can be taken as real- 
valued; under standard approximations, the tensor g is a 
quantity symmetric with respect to the exchange of any 
couple of indexes. 

By letting a s (t) = A s (t) exp [iip s (t)], the H can be 
rewritten as 



E 

{spqr} 



H(G, if) = H Q + V] G spqr cos(cp s + tp p -(p q - ip r ) (6) 



where H Q = J2 S V 8 (A S ) only depends on the amplitudes 
and G K „ar = 2gs P qrA s A p A q A r . As discussed in the liter- 
ature [n( Eqs. ^ are Langevin equations for a system of 
iV particles moving in 2N dimensions and the invariant 
measure is given by exp(— TL/kBTbath)- 

In a standard laser the resonant frequencies are equis- 
paced and this gives rise to various formulations of laser 
thermodynamics, which are based on the fact that the 
lu s + lu p = uj q + uj r will o nly involve a limited number 
of interacting modes [T3, El- The situation is drasti- 
cally different for nano-structured systems displaying a 
photonic-band gap. It is indeed well established that in 
proximity of the band-edge, a DOS enhancement with 
respect to vacuum is obtained. All the corresponding 
modes will have overlapping resonance such that lo s = loq, 
(where ujq is the position of the peak in the density of 
states, which is assumed to be in correspondence of the 
resonance of the amplifying atomic medium); addition- 
ally the resonance condition lu s + uj p = uj q + w r need not 
to be exactly satisfied, but it sufficient that linewidth of 
the corresponding modes needs to be overlapped for a 
relevant interaction (2?]]- Hence for such a kind of sys- 
tem it is interesting to consider a mean field regime where 
all the modes interact in a limited spectral region around 
u)q. For the mode- locking transition one can limit to con- 
sider the phase dynamics. Indeed ML entails the passage 
from a regime in which the mode-phases are independent 
and rapidly varying ( "free-run" regime or "paramagnetic 
phase" in the following) [25] on times scales of the order 
of 10 fs 28] , to a regime in which they are all locked at the 
same values ( "ferromagnetic phase" ) . In correspondence 
of this transition the laser output switches from a contin- 
uous wave noisy emission to an highly modulated signal 
(which is a regular train of short-pulses for equi-spaced 
resonances). Mode-amplitude dynamics is not affected 
(at the first approximation) by the onset of ML. Indeed 
for lasing modes 7 S > a s , so that the potential V S (A S ) 
has a single minimum in A s = \J {pf s — a s )/g s . Thus the 
amplitudes of the lasing modes will fluctuate around this 
minimum: we neglect these fluctuations, which are small 
if the potential well is deep enough, and threat the A s 
as quenched variables (£ s — uj s \a s \ 2 = (£ s ) = loqA 2 for 
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the relevant modes). The relevant subspace spanned by 
the system is given by the phases, that are taken as the 
dynamic variables (see, for example, [l[ for a discussion 
of the role of mode-phases with respect to amplitudes in 
ML processes). 



III. LOCALIZED VERSUS DELOCALIZED 
MODES IN THE SELF-MODE-LOCKING 
TRANSITION 

In previous works [Lsl . [l^ | we made reference to a com- 
pletely random resonator, for which the G coefficients 
were taken Gaussian distributed with zero mean. This 
is the natural approach when dealing with strongly dis- 
ordered resonators, in which the involved modes can be 
localized or delocalizcd in the structure and the corre- 
sponding resonances and spatial distribution can have 
different degrees of overlaps (as in [2!|). Here we make 
reference to the opposite regime, corresponding to case 
in which the structure is quasi-ordered, with the pres- 
ence of a small amount of disorder. The disorder is such 
that the variations in the coupling coefficients g sp q r can 
be taken as negligible with respect to their statistical av- 
erage (g S pqr) — g, however it is sufficient to induce the 
existence of a tail of localized modes in the photonic band 
gap [U. 

As discussed above, we consider mode-resonances 
packed in a small spectral region Alu ( if compared with 
the central carrier angular region, i.e. Alu << luq). This 
kind of system is very different from the standard laser 
cavity, with equispaced mode- frequencies. A prototypical 
structure is given by a photonic crystal doped by active 
materials. In the absence of disorder the involved modes 
are Bloch modes, which are extended over the whole sam- 
ple and are absent in the forbidden band gap. Their DOS 
is peaked at the band-gap edge ||. In this case, the 
modes have overlapping resonances (the width of each 
spectral line being determined by material and radiation 
losses) and they also have not-negligible spatial overlap 
(with exception of those mode-combinations which are 
vanishing for symmetry reasons). In the presence of a 
small amount of disorder it is well established [3(| that 
a tail of localized states appears in the photonic band 
gap. Hence the localized states also have overlapping 
resonances in tiny spectral regions in proximity of the 
band-edge. Their spatial overlap can be strongly reduced 
with respect to the Bloch modes, however the exponen- 
tial tails of their spatial profiles are expected to provide 
not vanishing values for g. Localized states can also be 
introduced intentionally, e.g. by using defects in a planar 
PC slab- waveguide [1], or in coupled cavity systems (see 
e.g. [f| and references therein). In this case the spatial 
overlap and the resonance frequencies can be tailored at 
will. Thus, in the general case, extended and localized 
states can be involved in the mode-locking transitions 
here considered. However, simple scaling arguments lead 
to the conclusion that the two kind of modes display a 



macroscopic difference in their "thermodynamics" . 

Extended modes. We start considering the extended 
modes. From the normalization, the modules of the 
eigenvectors E s are such that J dV\E s \ 2 — const. De- 
noting with V the volume over which a given mode is 

— 1/2 

different from zero, one has E s ~ V . For extended 
modes V oc V; in addition, most of the mode-overlaps 
are not vanishing, hence the sum in ([6]) will involve all the 
modes (within the spectral region Aw), and this implies 
that H] ex V 4 oc N 4 . The Hamiltonian is hence more 
than extensive (for an extensive one H oc N); however a 
thermodynamics approach is still possible if one accepts 
that the effective temperature will depend on the volume: 
the effective temperature will be taken as proportional to 
TV -2 oc V~ 2 , as detailed below, so that the resulting ef- 
fective Hamiltonian will be proportional to V. Physically 
this corresponds to the fact that the energy (and hence 
the pumping rate) needed to induced the ML transition 
will grow with the number of modes, if only extended 
modes are involved. For the overlap coefficients, in the 
absence of strong disorder g spq r —9 7^ 0, and the cou- 
pling g scales with the inverse volume, g oc V~ x oc iV _1 
(\E\ 4 ~ V~ 2 and J dV\ ~ V). Thus considering the in- 
variant measure (^4 S ~ A, g sp qr — g 6 G spqr ~ G = gA ), 
one has 



exp(-- 



ksTbgth 

eX P[~T^F C0S (^ + Vp ~ - fr)] = 

KB 1 bath ^ 
spqr 

. gA 4 N 3 1 „ 

eX Ph-r^ — 2^ C0S (^s +<Pp- <Pq~ <Pr)\ 



(7) 



exp[ 



KBTbath N 3 



where (3 = \g\A 4 N 3 /k B T bath = 1/T oc iV 2 is an inverse 
adimensional temperature, and the mode-phase depen- 
dent (extensive) Hamiltonian is given by (within an ir- 
relevant additive term) 



H (ext) = £ [i - cos (^ 

spqr 



+ tp p -(p q - tpr)] . (8) 



In ([8]) we have used the fact that g < in all the phys- 
ically relevant regimes. No transition is expected for 
g > 0, as detailed below. Note that, as we assumed that 



the condition s +0 t 



j+0 r is not exactly satisfied, and 



possibly due to the presence of disorder, the integrand in 
might have oscillations or fluctuations in sign that can 
affect the scaling of g with volume (the case of completely 
random fluctuations lead to a scalin g g oc V~ 3 / 2 and to 
an extensive Hamiltonian as in [l8l. Il9l|). Intermediate 
regimes might be present depending on the strength of 
the disorder. 

Localized modes. Next we consider localized modes, 
that exponentially decay in space. Defining V as above, 
it turns that it does not scale with the volume V of the 
sample, but can be written as V a = L 3 where L is an av- 
erage localization length. The overlap coefficients g avqr 
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will have a statistics strongly peaked around some av- 
erage value g, measuring the average amount of spatial 
overlap between localized modes that are neighborhood 
in space. Hence the sum in the Hamiltonian will only 
involve first neighborhoods and Hi oc V, while g will not 
depend on the size of the system. For the invariant mea- 
sure we will have (the angular bracket in (spqr) denoting 
sum over first neighborhoods): 



where H is 



exp(- 



exp[ 



kBTbnth 

gA* 



) 



knT t 



B-Lbath 



exp[-/3H(' oc )] , 



[spqr) 



C0S((/? S + (p p - lp q - (p r )] 



(9) 



where (3 = \g\A 4 /ksTbath = 1/T will not depend on the 
system size (we stress that the system size must be such 
that a large number of localized modes are present), and 
the Hamiltonian (within irrelevant additive constants) is 

H (loc) = [1 - COs(p s + <p p - <f q - tp r )} . (10) 

(spqr) 

Summarizing, if the transition involves extended modes, 
the effective temperature for the critical transition is ex- 
pected to depend on the size of the system; conversely for 
localized modes the critical temperature will be indepen- 
dent on the system size. 

Unfortunately, at variance with fully-connected (or 
"mean field") models as ©, analytical treatment of 
short-range Hamiltonians as (jTUJ) is almost always im- 
possible and the analysis can only be numerically per- 
formed. However, it is well established within the sta- 
tistical physics community that mean-field models ob- 
tained from first-neighborhood systems conserve most of 
the thermodynamics properties, and more specifically the 
existence of a thermodynamic transition, at least above 
the so called "lower critical dimension" di (for example, 
for the Ising model is di = 1, for the XY model is di = 2). 
Our model falls in the class of XY models so we expect 
that the transition exists, as long as d > 2, also in the case 
of localized modes, and the following analysis applies at 
least qualitatively. Since ([5]) can be analytically treated 
for thermodynamic, topological and dynamic properties 
we will limit to this model in the following. The existence 
of a thermodynamic/topological transition is expected in 
the general case, while the different scaling properties of 
the effective temperature enables to discern localized and 
delocalized interactions. 



IV. THERMODYNAMICS 

In this Section we study the thermodynamics of the 
laser Hamiltonian in the mean field approximation ([5]), 
within the quenched amplitudes approximation. The 
partition function is 



Z = dip e 



(11) 



H 



-y 

spqr 



[1 - COs(ip s +(p p -(p q - p r )} . (12) 



The Hamiltonian (fT2"|) is very similar to that defining 
the fc-trigonometric model (fc-TM) for k = 4, introduced 
in [23j with the aim of studying the relation between 
phase transitions and topological property of the poten- 
tial energy surface [2l|, [3l|, [32| . 
Defining the "magnetization" 



(13) 



where £ and ip depend on {pi}, we have 



H = Re 



Re 



i ^ t 1 _ expi(p s + ip p -p q - p r )} 



spqr 



2 ) 



N(l-£ 4 ) 



(14) 



By definition a vanishing z denotes uncorrelated phase, 
as in the "free run" regime, conversely if z ^ the 
phase of the modes are correlated and locked. The ther- 
modynamics of the mean-field model is exactly solved 
by neglecting the correlations between different degrees 
of freedom, and obtaining an effective Hamiltonian that 
contains a parameter to be determined self-consistently. 
Introducing the mean (complex) "magnetization" £ = 
(e llp ), and substituting in equation (fl"2"]) the expression 



J(<Ps+<Pp-<P< l —<Pr) 



+ { e 3V.)( e 3<Pp) e -3V, ( e -JV-) + 



-3 (e J ^ ) (e^* ) (e-^« ) (e~^ ) = 

= 2CC* (e jip C + e~ Jlp C) ~ 3C 2 C* 2 = 4C 3 cos p - 3C 4 (15) 

where the last equality stands because we have chosen 
C to be real without loss of generality (corresponding to 
choosing a particular magnetization of the low temper- 
ature state), the effective Hamiltonian h per degree of 
freedom reads as 

h{p) = l + 3C 4 -4C 3 cos</? . (16) 

The self-consistent equation for £ turns out to be 



Ii(4/?C 3 ) 

c = {coaip)h = WfW) 



(17) 



where Iq(o) — (27r) _1 f Q * dp) exp(a cos tp) and I\(a) — 
Jg(a) are the modified Bessel function of order and 1, 
and (. . . )h is the average over the probability distribution 



P(p) 



-0Hv) 



(18) 
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FIG. 1: (Color online) Free energy /(C) - /(C= 0) as a function 
of magnetization £ at different temperatures. From high to 
low: T = 0.910, T = 0.717, T = 0.616, T c = 0.548, T = 
0.504. T marks the appearance of the unstable minimum, T c 
is the transition temperature at which the solution with £ > 
becomes thermodynamically stable. 



with 



d(p e -f3h{ 9 ) = 2n e -/3(i+3C 4 ) / ( 4/ 3C 3 ) . (19) 



The solution of Eq. (fT7|) are the extrema of the free 
energy / as a function of £ 



0f = -lnZ = (3(1 + 3C 4 ) - In2^ (4/3C 3 ) 



(20) 



whose absolute minimum is the thermodynamical stable 
solution. 

The value £ = 0, corresponding to the paramagnetic 
solution, always solves equation (jTTJ) , but it gives the 
stable (lower free-energy) solution only for low j3 (high 
T). On lowering T, at T Q = 0.717 other solutions appear, 
such that C 7^ 0- However, the stable solution is still the 
paramagnetic one £ = 0. At T c — 0.548 the solution 
( 7^ becomes the stable one, and a first-order phase 
transition takes place. In Fig. [T]the ^-dependence of the 
free energy / is reported for different temperatures. The 
stable solution C(T) is shown in Fig. [5^ (full line) while 
dashed lines denote unstable solutions (local minimum 
and maximum). In Fig. [2}d the T-dependence of the 
energy 



(21) 



is shown for the stable (full line) and unstable (dashed 
lines) solutions. 

A remark on the sign of the coupling g. For g > 
the sign of the cos term in the Hamiltonian (jT3J) 
is positive and the self-consistent equation reads £ = 
-Ji(/3A4C 3 )/7 (/3A4C 3 ), which has the only solution ( = 
0. Then, in this case the phase transition does not take 
place. 
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FIG. 2: (Color online) Magnetization £ - part (a) - and en- 
ergy e - part (b) - as a function of temperature. Full lines 
correspond to the thermodynamically stable solution of Eq. 
(|17[) . dashed lines are the unstable solutions. At T c — 0.548 a 
first-order thermodynamic phase transition takes place, while 
at T = 0.717 unstable solutions appear. 



V. TOPOLOGY 

After having ascertained the existence of a first-order 
phase-transition, we consider the property of the station- 
ary points (saddles) of the potential energy landscape 
of the system 1231. As said in the Introduction, in re- 
cent works [2l[ |22|, it has been conjectured that phase 
transitions are signaled by discontinuities in the config- 
uration space topology. More precisely, for a system de- 
fined by a continuous potential energy function V(q) (q 
denotes the N-dimensional vector of the generalized co- 
ordinates) a thermodynamic phase transition occurring 
at T c (corresponding to energy V c ) is the manifestation 
of a topological discontinuity taking place at V c ("topo- 
logical hypothesis"). The most striking consequence of 
this hypothesis is that the signature of a phase transi- 
tion is present in the topology of the configuration space 
independently on the statistical measure defined on it. 
Through Morse theory, topological changes are related to 
the presence of stationary points of V , and, more specifi- 
cally, to the discontinuous behavior of invariant quantity 
defined on them, as the Euler characteristic x- Subse- 
quent works [33l . l34j | have shown that, at least for some 
model system, a "weak topological hypothesis" applies in 
place of the "strong" one: the Vg at which a topological 
transition takes place does not coincide with the thermo- 
dynamic one V c 7^ V$ , but is related to it by a saddle-map 
M, from equilibrium energy level to stationary point en- 
ergy: M(V C ) — Vg. Then, the role of saddles has been 
demonstrated to be of high relevance for the topological 
interpretation of thermodynamic transitions. Here we re- 
port on the saddle properties of the considered nano-laser 



model. 

The stationary point Tp are defined by the condition 
&H(Tp) — and their order is defined as the number 
of negative eigenvalues of the Hessian matrix fly = 
(d 2 H / d(fid(pj)\ip. To determine the location of station- 
ary point we have to solve the system 



OH 



4^ 3 sin(^-V) = ) 



(22) 



where we have used Eqs. (fT3j) and (fT4|) . 
A first group of solutions arises for £ = 0; from equa- 
tion flUl) we have H = N(l - £ 4 ), and then the sta- 
tionary points with £(Tp) = are located at the energy 
e = H(Tp)/N = 1. Now we restrict ourselves to the region 
e ^ 1 because, as we will see at the end, the quantities 
in which we are interested are singular when e = 1. For 
e 7^ 1, equation (j22|) becomes 



sin^fc - tp) = 
and its solutions are 

Pk = [tp + m k Tr] 



mod 2tt 



(23) 



(24) 



where irik = {0,1}. The unknown constant if> is found 
by substituting Eq. in the self-consistency equation 



Introducing the quantity n(Tp) defined by 

n = N- 1 ^2m i , l-2n = iV- 1 ^(-l) r 



we have, from equation (|25p . 

£ = 1 - 2n . 



(25) 



(26) 



(27) 



As £ is positive defined, the only solutions are for n < 
1/2: there are not stationary points with n > 1/2. For 
n < 1/2 tp can assume all the values in [0, 2tt) for any 
choice of the set {m^} and all the stationary points of 
energy e ^ 1 have the form 

WT = bl> + "^Lod 2x - where m = i m k}, m k = (0, 1) 

(28) 

under the condition 



n = N^^mk < 1/2 . 

k 

The Hessian matrix is given by 

Hij = 4£ 2 { - -^f 2 sin(<^ - ip) sm(<^- - ?/>)- 
^■cos(<^i — <Pj) + 8ij£cos(tpi - ip)\ . 

In the thermodynamic limit it becomes diagonal, 



Hi, 



■i^Sij COs(( y 9. i - -0) 



(29) 



(30) 




FIG. 3: (Color online) The entropy of saddles a as a function 
of potential energy e. 

Neglecting the off-diagonal contributions (their contribu- 
tion changes the sign of at most one of the N eigenval- 
ues (23[) the eigenvalues of the Hessian calculated at 
the stationary point Tp are obtained substituting equation 
(El in ffl), 



A fe = (-1)™*4£ 3 



(31) 



Therefore, the stationary point order v(jp), defined as the 
number of negative eigenvalues of the Hessian matrix, is 
simply the number of mk = 1 in the set m associated with 
Tp; we can identify the quantity n(Tp) given by equation 
(gj|J) with the fractional order v(Tp)/N < 1/2 of Tp, Then, 
from equation ((14)) and (f27|) we get a relation between 
the fractional order n(Tp) and the potential energy e(Tp) = 
H(Tp)/N at each stationary point Tp. It reads 

1 r 



1 



1 



a/4 



(32) 



where we have used the condition n < 1/2. 
Equation (f3"2"|) brings the condition 

1 - e > 0, (33) 

so there are no stationary points for e > 1, while for e < 1 
the fractional order n — v/N of the stationary points is a 
well defined monotonic function of their potential energy 
e, given by equation (|32|) . 

The number of stationary points of a given order (apart a 
degeneracy factor) is proportional to the number of ways 
in which one can choose v times 1 among the {mk}, i.e. 
(1). Following [H, its logarithm 



a{e) = liniTv^oo j? In 
— n(e) In n(e) 



N \ _ 
Nn(e)) 

(1 - n(e)) lnfl - n(e)) 



(34) 



represents the configurational entropy of the saddles. 
Substituting in this expression equation (|32p we have 



1 



<e)=-- l-(l-e) 



a/4 





'1 " 


In 






2 . 



l-(l-e) 



i/4 
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1 + a-e) 1 / 4 





'1 " 


In 






2 . 



l + (l-e) 



1/4 



(35) 



For e > 1 indeed we have, obviously, <j(e) = 0. This 
quantity is related to the Euler characteristic x °f the 
manifolds M e = {ip\H(ip) < Ne} [23| and its singular 
behavior around the point e = 1 is related to both the 
presence and the order of the phase transition that oc- 
curs. 

In figure [3] the quantity a is reported as a function of 
energy e: one can see that the presence of a phase transi- 
tion is signaled by a singularity of o~(e) at the transition 
point e = 1 . It is worth noting that the curvature of the 
quantity <r(e) around the transition point e — 1 is posi- 
tive, according to what found in Ref. [23| for first order 
transitions. 

Summarizing, the study of stationary points shows 
that the presence of the phase transition is signaled by 
the stationary point properties. More specifically, the 
singular behavior of the configurational entropy cr(e) at 
transition point e = 1 is the topological counterpart of 
the thermodynamic transition. We note that these find- 
ings do not allow to discriminate between the strong and 
weak topological hypothesis, as in this case the map M 
from equilibrium energy levels to stationary point ener- 
gies is triv ially the identity at the transition point e = 1 : 
M(l) = 1 [Hf. 



VI. COHERENCE PROPERTIES AND 
DYNAMICS 



The dynamics of interacting lasing modes close to the 
mode-locking transition can also be investigated; it leads 
to explicit results for measurable correlation functions. 



A. Single-mode first order coherence 

We start considering the single- mode ( "self" ) first or- 
der coherence : 

5 " W - (E*(to)E n (t )) ' (36) 

where E n (t) — ^/uJ^a n exp(— iu n t) is the electric field 
emitted at the angular frequency tv n (omitting an 
inessential factor depending on the point in space where 
the field is measured). For quenched amplitudes, = 
luoA 2 , one has (omitting the mode- index n) 



(l) (t) = F {1) (t) -iupt 

9 [) FW(0) 



(37) 



where (t) is the unnormalized single mode first order 
coherence given by 



and we used the fact that at equilibrium the time aver- 
age over t can be replaced by a statistical average [35j . 
Below the threshold (T > T c ) the phases are uniformly 
distribute in [0, 2tt), so that F^'(t — > oo) = 0. In con- 
trast, at mode-locking the phase is blocked around a fixed 
value, thus F^(t -> oo) ^ 0. 

We are interested in the time-delay-profile of coherence 
function given by 



F(t) = FW(i) -F^{oo) 



(39) 



F(t) can be explicitly calculated in the mean field theory 
(see Ref. [24[ for all the details of the computation) . In 
fact it can be shown that the single mode dynamics can 
be mapped into the effective equation 



jip(t) = -4C 3 smip(t) + r)(t) 



(40) 



where £ is the thermodynamic value of the magnetization 
we determined in section (|IV[) . -q a ^-correlated Gaussian 
noise with variance 27T, and 7 a constant fixing the time- 
scale (in the following we use units such that 7 = 1). 

From Eq. (j4"0|) it is evident that in the paramagnetic 
phase (( = 0) the phases will freely diffuse, while in 
the ordered (mode-locked) phase they fluctuate around a 
given value that, without loss of generality, can be taken 
as (p = 0. Eq. (|4"0"]) can be solved and the self-correlation 
function of a single mode 



F(t) 



( e Mt) e - 



i<p(0)\ 



J<p(t) 



He 



> (41) 



can be computed [24j |. Using symmetry properties, the 
above function can be written as 



F(t) = F c {t) + F s (t) , 



where 

Fc(t) 
F s (t) 



(cos (fi(t) cos </?(0)) 
(sin ip(t) sin</j(0)) . 



(cos (p(t)} (cos 1^(0)) 



(42) 



(43) 



(38) 



Following [24] the self correlations (|43|) can be numeri- 
cally determined and they turn out to be nearly expo- 
nential at all temperatures, F c (t) oc e _t / Tc and F s (t) oc 
e -t/r._ j n piggj the function F(t)/F(0) is shown for dif- 
ferent temperatures. Upon increasing T the decorrelation 
time increases for T < T c , while it decreases for T > T c 
(after a sudden jump at T c ). This behavior is evident an- 
alyzing the T-dependence of the relaxation time. In the 
upper panel of FigJS] the quantity t c is shown as a func- 
tion of temperature. Full lines refer to stable states, while 
dashed lines to unstable ones. We note that in paramag- 
netic high-T phase r c (and r s ) have a 1/T-dependence, as 
expected for free Brownian motion. In the low-T phase 
the behavior of r s (not showed in the Figure) is very sim- 
ilar to that of t c . 

Summarizing, in absence of mode-locking the single- 
mode first-order coherence function has an exponential 
trend (corresponding to a Lorentzian linewidth), whose 
relaxation time decreases as the average energy per mode 
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FIG. 4: (Color online) Time dependence of the self correlation 
function F(t) at different temperatures. From left to right: 
T = 0.1,0.5,0.6,1.0 (T c = 0.548). Time is in normalized 
units. 



is reduced (i.e. the temperature is increased). At the 
mode-locking transition, the coherence function is ex- 
pressed as the sum of two exponentials (corresponding 
to the two quadratures of the phase-modulated laser sig- 
nal), whose time-constants have a jump with respect to 
the "free-run" regime, and decreases while increasing the 
average energy per mode (and hence reducing the tem- 
perature). 



B. Multi-mode first order coherence 

Here we consider the multi-mode ("collective") first 
order coherence: 



(E*(t )E(t +t)} 



(E*(t )E(t Q )) ' 



(44) 



with E(t) = ^ n y/Lu^a n exp(—iui n t). In the quenched 
amplitudes approximation, proceeding as above, it is pos- 
sible to write 



G«(0) 



(45) 



where we have taken u n = loq (since all the modes are 
taken as densely packed around loq, small differences be- 
tween u) n and wo can be embedded in the phase tp n ) , and 

G(1) W = ]^£< eW * )e " W0) > = (*(*)**(0)) . (46) 

nm 

is the correlation function of the magnetization z(i) = 
N^ 1 J2 n ex P IWn(*)]- We can write 

G«(t) = GW(oo) + iG(t) = C 2 + ±G(t) , (47) 

where C?W(oo) = C 2 is the asymptotic value (we recall 
that £ = (z) is assumed real), which is acquired at the 
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FIG. 5: (Color online) (a) Relaxation time r c of the self corre- 
lation function F c (t) as a function of T. Full line corresponds 
to stable state (paramagnetic above T c and ferromagnetic be- 
low it). Dashed lines refer to unstable solutions, (b) Relax- 
ation time r^f of the collective correlation function G c (t) as a 
function of T. Full and dashed lines as before. 



mode-locking transition, and the collective connected cor- 
relation function G(t) is defined as 



G(t) = N[G^(t)-G^(oo)] = 
1 



nm 



itp n (t) 



He 



-i<p m (0)\ 



(48) 



This function has a finite limit for N — ► oo, see the Ap- 
pendix, which can be computed following [24[. As for 
F(t), G(t) can be written as a sum of two terms 



with 



G(t) = G c (t)+G s (t) , (49) 
G c (t) = y^[(cos y> w (t)cosy m (0))- 

nm 

- (cos <p n (t)) (cos <p m (0))] , (50) 
G„(t) = -^^[(sin^„(i)sin^ m (0)) . 



Again one finds an exponential decay, G c (t) oc e { / r ° 
and G a (t) oc e _t / T = [24] . In the lower panel of Figf5]the 
relaxation time is shown as a function of tempera- 
ture. We hence expect a trend for G(t) which resembles 
F(t); however it is worth noting that the quantity t~ 
diverges when T — > T~ , that is, when, starting from 
low T phase, the point where the unstable solution 
(dashed line) disappear is approached. The occurrence 
of the thermodynamic transition at T c < T Q prevents the 
divergence of . The (not shown) does not diverge 
at T„. 
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Multi-mode second order coherence 



VII. CONCLUSIONS 



We also consider the multi-mode (collective) second- 
order coherence: 



P)m _ (E*(t )E*(t + t)E(t )E(t + 1)) 



again with 

E(t) 

We get 



(E*(t )E(t )) 2 
(I(t )I(t +t)) 
(I(to)) 2 ' 

= v / ZJ^"a n e _2a;Tli oc z[t)e 



,(2) 



(*) 



G( 2 )(oo) 



(51) 



(52) 



(53) 



where G^(£) = (2(i)z*(t)z(0).z*(0)). We can decompose 
this function in its connected components (see Ref . [36j| . 
Eq. 4.23); recalling that ( — (z) is assumed real we get 



G< 2 >(t) = (z(i)z*(t)z(0)z*(0)) c + CG 3 (t) 
+ (z(0)z*(0)) 2 + |<z(i)z(0)) c | 2 
+ |(z(t)z*(0)) c | 2 + C 4 , 



(54) 



where the function Gs(t) is a sum of connected 
three-point functions: G 3 (t) = 2Re[(z(t)z*(t)z(0)) c + 
(z(t)z(0)z*(0)) c ]. Using the results of the Appendix, we 
have (z(t)z*{t)z(0)z*(0)) c oc N~ 3 and G 3 oc N~ 2 ; from 
Eq.s (JSUJ), we have 



(55) 



(z{0)z*{Q)) c = N- l G{Q) , 
(z{t)z(0)) c = N-^G^t) - G s (t)} , 
(z(t)z*(0)) c = N-^Gdt) + G s (t)} . 
Then we obtain 

G(0) 2 + 2G c (t) 2 + 2G s (t) 2 



G< 2 )(i) = C 4 + 



N 2 



(56) 



+ (G 3 (t) + 0(N- 3 ) 



In the paramagnetic phase, £ = 0, and by symmetry 
G s (t) — G c (t) — G(i)/2; moreover these functions tend 
to for t — ► oo. Then we get 



9 {2 \t) 



G{Qf + 2G c {tf + 2G s {tf 



= 1 



G(oy 



G(0) 2 

■ = 1 + |5 (1) W| 2 , 



(57) 



where we make use also of Eq. (1471) . This result is indeed 
what we expect for light modes evolving independently 
and rapidly. 

In the mode-locked phase Eq. (|5"Tj) will not hold but 
a relation between and can still in principle 
be deduced from the knowledge of function G 3 (£), us- 
ing Eq. (H)6"l) and PT)) . Future works will address this 
point. 



By using a simple model that is expected to describe 
multi-mode dynamics of tightly packed extended and/or 
localized modes in a nano-optical resonator, we predict 
the existence of a first-order phase-locking transition 
when the averaged energy per mode is above a critical 
value (correspondingly the adimensional effective tem- 
perature is below T c ). This value depends on the aver- 
age value of the mode-overlap coefficient g. If the transi- 
tion involve extended modes, one has (omitting indexes) 
g = oj 2 J \E\ 4 dV = XoUqV -1 , with xo a reference suscep- 
tibility value. Conversely if localized modes are involved 
it is g = xoVd 2uj o wnere Vq is the average localized mode 
volume (Vo cx h\ with Lq the localization length). In 
the former case, for a fixed spontaneous emission noise 
Tbath, it is found that the critical mean energy per mode 
is U-dependent 



S^=u;oA 2 = J kBTbathV , 

V T eX0 



(58) 



while for localized modes 



£^ = uj A 2 = 



' ksTbathLQ 

TcXo 



(59) 



Hence the critical energy for the phase-locking transi- 
tion has very different scaling behavior with respect to 
the system size, depending on the degree of localization 
of the involved modes. In the general case one can ex- 
pect intermediate regimes between those considered, so 
that the trend of the critical energy (determined by the 
amount of energy pumped in the system per unit time, 
i.e. the pumping rate) versus system volume is an inter- 
esting quantity which can experimentally investigated. 
A topology-thermodynamics relationship has been evi- 
denced for our model, corroborating previous findings on 
this topic: the thermodynamic transition is signaled by 
a singularity in the topological quantity a. 

The exact solution of the dynamics of the model pre- 
dicts the divergence of the relaxation time of the first- 
order coherence function g"> at the transition. This be- 
havior might be observed in experiments; the different 
scaling with respect to the systems size also affect the 
position of the transition, as determined by the predicted 
jump in the relaxation time or, equivalently, in an abrupt 
change of the single-mode laser linewidth while varying 
the pumping rate. 

This analysis points out the rich phase-space structure 
displayed by these systems, while varying the amount of 
disorder or the profile of the density of states. Hence 
nano-lasers not only may furnish the basis for highly in- 
tegrated short-pulse generators, but are a valuable frame- 
work for fundamental physical studies. These deserve fu- 
ture theoretical and experimental investigations and can 
be extended to other nonlinear multi-mode interactions. 
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Appendix 

We discuss here the scaling of the correlations of z 
when TV — > oo. 

The basic fact is that the variable z is intensive and, for 
a mean field system, has a probability distribution of the 
form [H 



This simple argument shows that (z k ) c oc iV 1 k . In par- 
ticular, 

(z 2 ) c = (z 2 ) - (z) 2 oc AT- 1 , 

(z 3 ) c = (z 3 ) - 3(z)(z 2 ) c - (z) 3 oc N- 2 , (63) 
(z\ = (z 4 )-A(z)(z 3 ) c -3(z 2 ) 2 c -(z)^N- 3 . 



Pn{z) 



(60) 



where F is related to the thermodynamic free energy of 
the system. Consider the generating functional 



z NF{j) = {e N ]z) = / dze N[F(z)+jz] ^ 



(61) 



then for large TV, F(j) = ma,x z [F(z)+jz] and it is a quan- 
tity of order 1. The connected correlation functions [36[ 
are derivatives of NF(j) with respect to Nj: 



S k NF{j) 



(SNjY 



= N 



3=0 



(SjY 



(62) 



3=0 



We assumed that z is real but the same derivation can 
be repeated for a complex variable; the only difference is 
in the definition of the connected correlation functions. 

For the dynamics, one can write a similar expression 
for the probability of a trajectory z(t), see [24[ and ref- 
erences therein: 



P N [z(t)] = e NF Wn 



(64) 



Repeating the derivation above using functional integrals 
one obtains exactly the same results for the scaling with 
N. 
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